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Abstract: We generalize the Milne quantization condition to non-Hermitian systems. 
In the general case the underlying nonlinear Ermakov-Milne-Pinney equation needs to be 
replaced by a nonlinear integral differential equation. However, when the system is PT- 
symmetric or/and quasi/pseudo-Hermitian the equations simplify and one may employ 
the original energy integral to determine its quantization. We illustrate the working of 
the general framework with the Swanson model and two explicit examples for pairs of 
supersymmetric Hamiltonians. In one case both partner Hamiltonians are Hermitian 
and in the other a Hermitian Hamiltonian is paired by a Darboux transformation to a 
non-Hermitian one. 


1. Introduction 

As one of the first phase amplitude methods Milne provided in 1930 |j]] a relation between 
the time-independent Schrodinger equation and a non-linear integrable equation referred 
to these days as the Ermakov-Milne-Pinney (EMP) equation [^, [lL |] or variants thereof. 

Solving either of the two equations for any generic energy will provide a solution for the 
other. In addition, the interrelation involves an auxiliary equation whose solutions lead 
to the exact energy quantization in a very general fashion. It should be emphasized that 
the Milne quantization is exact and the more popular WKB-approximation is obtained as 

a limiting case when the second order derivative term in the EMP-equation is neglected. 
While the latter method has been generalized [Q] to non-Hermitian PT-symmetric systems, 
this task has not been carried out for the more general Milne quantization procedure. 
The main purpose of this manuscript is to perform the first step in this direction and to 
demonstrate that a successful application of the Milne quantization procedure is indeed 
possible. 
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We analyze two types of non-Hermitian systems, in one case we exploit the fact that 
the model is quasi/pseudo-Hermitian, the Swanson model, and in the other that it is 
PT-symmetric, a supersymmetric pair in which one of the partner Hamiltonians is non- 
Hermitian. 

Our manuscript is organized as follows: In section 2 we recall the key features of the 
Milne quantization procedure and generalize it to a general non-Hermitian setting. In 
section 3 we discuss the Swanson model and in section 4 we provide two explicit examples 
for pairs of supersymmetric Hamiltonians, where in one case both partner Hamiltonians 
are Hermitian and in the other only one of them. Our conclusions and outlook are stated 
in section 5. 

2. The Milne quantization for Hermitian and non-Hermitian systems 

We commence by briefly recalling the key idea of the solution procedure and quantization 
method proposed originally by Milne in 1930 [|lj. Its starting point is the time-independent 
Schrodinger equation in the form 


if"{x) + k 2 (x)ip(x) = 0, 


( 2 . 1 ) 


where the continuous energy parameter E and the potential V{x) are combined into the 
local wavevector k 2 (x) = h 2 /2m[E — V(x)]. Assuming the solution to equation (2H) to be 
of the general form 

if{x) = N p(x) sin [4>{x) + a], (2.2) 

with normalization constant N, constant phase a, amplitude p{x) and variable phase 4>{x) 
a direct substitution leads to the constraining equations 


p"(x) + k 2 (x)p{x) = 


A^ 


p 3 (x)’ 


and p 2 (x)<f)'(x) = A, 


(2.3) 


with A being some arbitrary constant. The first equation in is known as the Ermakov- 
Milne-Pinney (EMP) equation 111 From ( |2.2|) it is clear that its solution together with 
a solution for the auxiliary equation for the phase function will lead to an exact solution for 
the time-independent Schrodinger equation (|2.1|) for generic values of E. Notice that when 


we neglect p"(x), the two equations in (24S) combine into 4>'(x) = k(x) which corresponds 
to the WKB approximation. In what follows we will employ Pinney’s Q general solution 
for the EMP-equation 1 


p{x) = J if\(x) + T-V’K*). 


(2.4) 


with p(x o) = p 0 ^ 0, p\x o) = p' 0 , — oo < xo < oo. Here 4> 2 are the two fundamental 
solutions of the Schrodinger equation (2.1) and W := 2 ) = ~ denotes 


the corresponding Wronskian. Integrating the second equation in (2.3) directly and taking 


1 There exist other types of solutions, such as for instance the one reported in jH| involving two free 

constants, which we may, however, suitable chose. 
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the initial conditions to be -i/h( x o) = 1, 0 ) = 0, ipi(x 0 ) = 1, ^4(xo) = A implies W = A 

and leads to the general solution of equation ( |2.1[ ) expressed in terms of the solutions to 


the EMP-equation 


= Np(x) sin 


W 



s + a 


(2.5) 


Next we implement the boundary conditions. Demanding the wavefunction ip(x) to vanish 
at the boundaries then implies the quantization condition 


m 


W (E) 


7T 


f 


p (s, E)ds = n € N, 


( 2 . 6 ) 


when p(x) is non-vanishing, meaning that any solution E n to I(E n ) = n constitutes a 
bound state energy. Note that the value of 1(E) is not sensitive to the normalization 
factors in the fundamental solution. 

When the potential and possibly also the energy eigenvalues are complex the general 
treatment is more involved. In that case we can make the Ansatz 


xf(x) =Np(x)e i ^ x \ (2.7) 

with p(x),4>(x) € R and separate the wavevector into its real and imaginary part k 2 = K+ir. 
The substitution of (2/7) into the time-independent Schrodinger equation then yields the 
two constraining equations when reading off the real and imaginary parts 


p"(x) + k(x)p(x) = p(x)(f>' ( x), and f"(x)p(x) + 2 (x)p'(x) + r(x)p(x) = 0. 
Combining these two equations generalizes the EMP-equation ( |2.3|) to 

1 


p (x) + n(x)p(x) = 


p 3 (x) 


A — / T(s)p z (s)ds 


with 


4>(x) = A / p 2 (s)ds — I p 2 (t) ( / r(s)p z (s)ds ) dt. 


-2, 


( 2 . 8 ) 


(2.9) 


( 2 . 10 ) 


Evidently when r = 0 we recover (|2.3| ). These equations are difficult to solve, even in an 
approximate fashion. However, we may assume that the quantization condition (2.6) still 
holds when W(E) € R and Im[p 2 (s, E)\ is an odd function in s. We will demonstrate below 
that these properties can be attributed to the VE -symmetry of the models. 


3. A quasi-Hermitian model, the Swanson Hamiltonian 

Quasi/pseudo-Hermitian Hamiltonian systems constitute a large subclass of non-Hermitian 
systems III They are characterized by the fact that their non-Hermitian Hamiltonian 
H can be mapped to an isospectral Hermitian counterpart h by means of a similarity 
transformation h = r/Hr / -1 . The map i] is sometimes referred to as the Dyson map || 
and satisfies certain properties. A prime example for which this map and all other relevant 
quantities are known in its explicit analytic form is the Swanson model © 

Hg = oj + 1/2^ + aa 2 + (3 , u,a,/3€. R, (3-1) 
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with a = y/uj/2x+ i/y/2ujp, of = \Ju/2x — i/y/2ojp. Evidently H$ is only Hermitian when 
a = (3, but its isospectral Hermitian counterpart is known to be |]TT[ 

hs = ^P 2 + !fx 2 , (3.2) 

with 


—A(a + /3)+w=F(a + /3 — Aa >)yjl - 


(1 db A)^ 1 

The eigenvalue spectrum for both Hamiltonians is 


A £ [—1,1]- 


E n = ^n + a/w 2 - 4a/3, n £ N, 


(3.3) 


(3.4) 


and thus real for uj 2 > 4 a/3. The corresponding time-independent Schrodinger equations 
are exactly solvable for both Hamiltonians. The two fundamental solutions for the one 
corresponding to hs can be expressed in terms of parabolic cylinder functions, but in the 
current context it is more convenient to employ the solutions in terms of the closely related 
Whittaker functions 


Vh7) 



E 


1 

4 



Q(x) -|—— M e 

V x 2 v //x_/z_|_ 


1 

4 



(3.5) 


ll>o(x) = —^W E 1 
y/x 2 v ^+’-4 



0(x) + A=W E 
y/x 


1 

4 



(3.6) 


We neglect here normalization factors for the above mentioned reason. Unlike the solutions 
in terms of parabolic cylinder functions this choice guarantees that ^ 1,2 ( x ) £ R or ^ 1,2 ( x ) £ 
iR, such that p(x),W(E) £ R. Using these expressions we compute the energy integral 
I(E) in ( [2.6D and depict our results in figure jl]. 

The energy eigenvalues are located precisely at the expected values at points of inflec¬ 
tion of the function 1(E). 


4. Non-Hermitian models with supersymmetric Hermitian counterparts 


Now we study a model in which we exploit the PT-symmetry of the system. We consider 
a pair of supersymmetric quantum mechanical [fi~2| . 13, 14, 15] models described by the two 
Hamiltonians 


H± = L ± L t = + U 2 (x ) ± U'(x) = -^ + V ± (x), (4.1) 

ax z ax z 

involving the so-called superpotential U(x). It is easily verified that the solutions to the 
time-independent Schrodinger equations H±i/j± = Eil>± are related to each other by means 
of the two intertwining operators L± 


L± : = ± ii + u(x) ’ 




1 

71 


F . 


(4.2) 
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E 


Figure 1: Energy integrals 1(E) for the Swanson model, with I((2n — l)/4-\/2) = I(E n _i) = n £ N 
for w = 1/2, a = 1/8, $ = 1/4 I((2n - l)/2 v / 2) = I(E n _ i) = n £ N for w = 1, a = 1/2, /3 = 1/4 
and I((2n — l)ll/4\/2) = I(E n _i) = n £ N for w = 3/2, a = 1, /3 = 1/3 


Denoting now the two fundamental solutions to the Schrodinger equation by if and Xi 
Ioffe and Korsch [ 161 found that the corresponding Wronskians and solutions to the EMP- 
equations 

W± := W (if±,X±) , P± = + X± (4.3) 

are related to each other as 


W+ = W-, and Epi = (. L±p 


(4.4) 


The first identity follows from a direct substitution of the wavefunction in ( |4.2| ) into the 
defining relation for the Wronskian, the use of the Schrodinger equation and recalling that 
dW/dx = 0. The derivation of the second identity follows from a direct evaluation. We 
also add here for later use an intermediate relation from that computation 


Ep\ = u 2 p 2 _ + U(p 2 _y + (O 2 + (x'_) 2 • (4.5) 


Let us now select our super potentials to be of a very specific type, such that one of the 
partner Hamiltonians is Hermitian whereas the other one is not. Such a setting allows us 
to test our assertions from section 2. Bagchi and Roychoudhury [|lj]] provided a necessary 
condition for such type of pairs and noted that one may even construct solvable models in 
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this case. Separating the real and imaginary parts in the superpotentials in the form 

U(x) = a(x) + ib(x), with a(x), b(x) G R, a(x) = ln&(x), 

2 ax 

they observed that one obtains a real and a complex partner potential 

3b' 2 b"(x) 


V-{x) = 
V + {x) = 


4 b{x) 2 
b"{ x) 


2 b(x 
b' 2 


— b(x) 2 € R, 

2 - b(x) 2 + 2 ib'(x) £ 


(4.6) 

(4.7) 

(4.8) 


2 b{x) 46(x) 

In the following it will be important to utilize the effect of the parity operator V and 
time-reversal operator T on the various quantities involved. Our main requirement is that 
V+ becomes VT -symmetric; which is achieved as follows 

VT : a(x) —a{x),b{x) —y b(x); VT : U(x) —$■ —U(x),V±(x) V±(x). (4.9) 

In order to obtain real eigenvalues E € M, usually referred to as the spontaneously unbroken 
VT -symmetric regime, we also require the wavefunctions to be symmetric with regard to 
the anti-linear PT-operator (l8|, |l9[ 

VT : ^±(®),-> ^±(x),x±(x) -)• X±(x),W±(x) -> -W±(x),p±(x) ->• p±(x). (4.10) 

When assuming that £ M, it follows from ( |4.5| ) and the subsequent use of the second 

relation in (4.4) that 


Im (Ep 2 + ) = Im [(L + p_) 2 ] = ^ (bp 2 _) . 


(4.11) 


This implies that for real energies there will not be any contribution to the integral in 
( |2.6| ) from the imaginary part of the integrand \/p\ as it will be an odd function. The 
assumption if_,x~ € M also guarantees that W— € M and therefore by the first relation in 
(4.4) W+ € R, which are the requirements mentioned at the end of section 2. 

4.1 A Hermitian/Hermitian supersymmetric pair 

As an illustration for the working of the conventional Milne quantization for supersym¬ 
metric pairs we first consider a well studied exactly solvable in the mathematical physics 
literature, [20, 21, 22], the Poschl-Teller model [^3[. Taking the superpotential to be of the 
form 

U(x) = Atanx — Kcotx, k, A € R, 0 < x < tt/2. (4-12) 


equation (4.1) yields the pair of potentials 

V±(x) = A(A ± 1) sec 2 x + k(k ± 1) esc 2 x — (A + k) 2 , 


(4.13) 


with V-(x) being the standard Poschl-Teller potential. The fundamental solutions are well 
known. We have 


Yq (x) = sin K x cos A x -zF\ 
i/jf{x) = sin 1 '' 1 x cos A x 2 T 1 


k + A — E k ^ E 


1 ■ 2 

;/«+-; sm x 


(4.14) 


1 — k + A — E 1 — he X E 3 .0 . . 

-,-;-k;sui x , 4.15 

2 2 2 ’ v ’ 
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and 


= sin K+1 xcos A+1 X 2 F 1 
1 P 2 (x) = sin“ K x cos A+1 x 2 -F 1 


2 + k + A — E2-\-k-\-\-\-E 3 .2 ,. . 

-2-,-j- -,K + -\sm 2 x ,(4.16) 

1 — K + \ — El — K + X+E 1 , 2 fAA>7\ 

- 2 ->- 2 -;r K;sm x ’ ( 4 - 17 ) 


where 2 -E 1 denoted hypergeometric function and we abbreviated E := 
Solutions to the EMP-equation are simply obtained from (2.4) 


/9 ± (x) = ^/[^f(*)] 2 + [V>2 (*)] 2 > 

which allows us to compute the energy integrals ([2d]) to 


w ± (E) r/ 2 _ 2 

Jo 


I±{E) = —^4 —L I p± 2 (s, E)ds. 


y/ {n + A) 2 + E. 

(4.18) 

(4.19) 


Our numerical computations of (4.19) are depicted in figure 0. 



E 


Figure 2: Energy integrals I±(E) for a supersymmetric pair of Poschl-Teller potentials for coupling 
constants k = 2, A = 3, with 7_(0) = 7 + (24) = 1, 7_(24) = 7+(56) = 2, /_(56) = 7 + (96) = 3, 
7_(96) = 7+(144) = 4, 7_(144) = 7+(200) = 5, 7_(200) = 7+ (264) = 6, 7_(264) = 7+(336) = 7, 
7_(336) = 7+(416) = 8 and 7_(416) = 9. 

For the selected values of the coupling constant k = 2, A = 3 the solutions to I±(E±) = 
n + 1 yield Eq = 0, E+ = E~ +l = 4(n + l)(n + 6) for n = 0,1, 2,... This is of course 
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the well known quantization condition obtained from demanding that lim x _>.o ^fix) = 
lirrij.^^/ 2 = 0, achieved by setting the first entry of the hypergeometric function 2 F 1 

to —n with n = 0,1, 2,... 

4.2 A Hermitian/Non-Hermitian supersymmetric pair 

Next we consider a superpotential giving rise to a Hermitian potential paired with a non- 
Hermitian potential as proposed in [17]. We take the superpotential U(x) to be of the 
form 

' ' Ael, (4.20) 


1 Z 

U(x) = —- tanh x + -(1 — 2A) sechx, 


such that the real and imaginary parts are related as in (4.6). As expected, when evaluating 

(4.21) 


(4.1) one of the partner potentials turn out to be real 

1 


V-(x) = -1 + (A — A 2 ) sech 2 x, 


whereas the other one becomes complex 

V+(x) = - — (l — A + A 2 ) sech 2 x + i( 2A — 1) sech x tanh.: 
albeit VT -symmetric. The fundamental solutions are in this case 


Vq (x) = sinh x cosh A x 2 -F 1 
iff{x) = cosh A X2T 1 i 


-; — sinh 2 x 


1- 


> 1 + 


1 1 

2 ’ 2 


; — sinh" x 


(4.22) 

(4.23) 

(4.24) 


and according to 


we obtain the solutions for the partner Hamiltonian as 


■*l>i( x ) = 


^t{x) = 


cosh A i (x) 

12 VE 


6 [2 cosh 2 x + (2A — 1) sinh x(sinh x — i)] 2-^1 


—; — sinh 2 x 


- - sinh 2 (2x) [4 E + 4A(A + 2) + 3] 2^1 


5 2 

//_ + 1, n + + 1; - sinh x 


(4.25) 


cosh A : (x) 


2(2A — l)(sinhx — i) 2 T 1 


1 


1 1 


~ 2 ’ 2 ;_sinh x 


1 3 


+ -,!+ + -; -; — sinh x 


(4.26) 


4 VE 

+ (l — 4 E — 4A 2 ) sinh x cosh 2 x 2 T) 
where fi . t := (2 + 2A ± \/l — 4B)/4. 

We have now all the ingredients to evaluate the energy integrals in ( 4.1£ ). Our results 
are depicted in figure [3]. 

For the selected values of the coupling constant A the solutions to I±(E^) = n +1 yield 
Eq = —42, E+ = E~ +l = — (n — 6)(n — 7) for n = 0,1, 2,... This is again the quantization 
condition obtained from demanding that lim a ;_^-|- 00 i/jf(x) = 0, achieved by setting the first 
entry of the hypergeometric function 2 T 1 to — n with n = 0,1, 2,.. .The remarkable feature 
is here that we can still use the standard formula for the Milne quantization even though 
one of the Hamiltonians is non-Hermitian. Notice that this feature can be attributed 
entirely to the VT -symmetry of the system, which is responsible for the vanishing of the 
imaginary part in the energy integral. 
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Figure 3: Energy integrals I±{E ) for a supersymmetric pair potentials V± in ( 4.21 ), ( |4.22| ) for the 
coupling constant A = 15/2, with /_(—42) = 7 + (—30) = 1, /_(—30) = /+(—20) = 2, /_(—20) = 
I+(~ 12) = 3, I- (—12) = /+(—6) = 4, /_(—6) = 7+(-2) = 5, J_(-2) = J+(0) = 6, and 7_(0) = 7. 


5. Conclusion 

We demonstrated that the Milne quantization procedure can be successfully adopted to 
non-Hermitian systems that are either quasi/pseudo-Hermitian or VT -symmetric. For each 
scenario we provided an explicit example. We proposed some generalized formulae for the 
generic non-Hermitian case, which are left as a challenge to be solved for some concrete 
example. 

Building on the success, it is to be expected that this method can be applied also 
to systems for which the quantization is still incompletely understood [^4|], such as the 
complex Mathieu system currently of great interest as it corresponds to the eigenvalue 
equation of the collision operator in a two-dimensional Lorentz gas. 
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